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1 Introduction 

The purpose of this report is to present the theoretical foundations of the inter- 
action of evanescent fields on an object. Evanescent electromagnetic waves are 
inhomogeneous components of the near field, bound to the surface of the scat- 
tering object. These modes travel along the illuminated sample surface and ex- 
ponentially decrease outside it [Born 99a, Nieto-Vesperinas 91, Mandel 95], e.g., 
either in the form of lateral waves [Tamir 72a, Tamir 72b] created by total inter- 
nal reflection (TIR) at dielectric flat interfaces, whispering-gallery modes in di- 
electric tips and particles [Hill 88, Owen 81, Benincasa 87, Collot 93, Knight 95, 
Weiss 95, Nieto-Vesperinas 96] or of plasmon polaritons [Raether 88] in metal- 
lic corrugated interfaces (see Section 5). The force exerted by these evanescent 
waves on particles near the surface is of interest for several reasons. On the one 
hand, evanescent waves convey high resolution of the scattered field signal, be- 
yond the half wavelength limit. This is the essence of near-field scanning optical 
microscopy, abbreviated usually as NSOM [Pohl 93, Paesler 96]. These fields may 
present large concentrations and intensity enhancements in subwavelength regions 
near tips thus giving rise to large gradients that produce enhanced trapping forces 
that may enable one to handle particles within nanometric distances [Novotny 97]. 
In addition, the large contribution of evanescent waves to the near field is the ba- 
sis of the high resolution of signals obtained by transducing the force due to these 
waves on particles over surfaces when such particles are used as probes. On the 
other hand, evanescent waves have been used both to control the position of a 
particle above a surface and to estimate the interaction (colloidal force) between 
such a particle and the surface (see Chapter 6) [Sasaki 97, Clapp 99, Dogariu 00]. 

The first experimental observation, demonstrating the mechanical action of 
a single evanescent wave (i.e., of the lateral wave produced by total internal 
reflection at a dielectric (saphire-water interface) on microspheres immersed in 
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water over a dielectric surface) was made in [Kawata 92]. Further experiments 
either over waveguides [Kawata 96] or attaching the particle to the cantilever of 
an atomic force microscope (AFM) [Vilfan 98] aimed at estimating the magnitude 
of this force. 

The scattering of an evanescent electromagnetic wave by a dielectric sphere 
has been investigated by several authors using Mie's scattering theory (addressing 
scattering cross sections [Chew 79] and electromagnetic forces [Almaas 95]), as 
well as using ray optics [Prieve 93, Walz 99]. In particular, [Walz 99] made a 
comparison with [Almaas 95] . Although no direct evaluation of either theoretical 
work with any experimental result has been carried out yet, likely due to the 
yet lack of accurate well characterized and controlled experimental estimations 
of these TIR observed forces. In fact, to get an idea of the difficulties of getting 
accurate experimental data, we should consider the fluctuations of the particle 
position in its liquid environment, due to both Brownian movement and drift 
microcurrents, as well as the obliteration produced by the existence of the friction 
and van der Waals forces between particle and surface [Vilfan 98, Almaas 95]. 
This has led so far to discrepancies between experiment and theory. 

In the next section we shall address the effect of these forces on particles from 
the point of view of the dipolar approximation, which is of considerable inter- 
pretative value to understand the contribution of horizontal and vertical forces. 
Then we shall show how the multiple scattering of waves between the surface and 
the particle introduces important modifications of the above mentioned forces, 
both for larger particles and when they are very close to substrates. Further, we 
shall investigate the interplay of these forces when there exists slight corrugation 
in the surface profile. Then the contribution of evanescent waves created un- 
der total internal reflection still being important, shares its effects with radiative 
propagating components that will exert scattering repulsive forces. Even so, the 
particle can be used in these cases as a scanning probe that transduces this force 
in a photonic force microscopy operation. 

2 Force on a Small Particle: 
The Dipolar Approximation 

Small polarizable particles, namely, those with radius a <C A, in the presence of 
an electromagnetic field experience a Lorentz force [Gordon 73]: 

F=(p-V)S + - c ^xB. (1) 

— * — * 

In Equation (1) p is the induced dipole moment density of the particle, and S, B 
are the electric and magnetic vectors, respectively. 
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At optical frequencies, used in most experiments, the observed magnitude 
of the electromagnetic force is the time-averaged value. Let the electromag- 
netic field be time-harmonic, so that S(r,t) = 3?e{E(r) exp(— tut)}, B(r,t) = 
Ke{B(r)exp(-iu;t)}, p(r,t) = 3?e{p(r) exp(-iut)}; E(r), B(r) and p(r) be- 
ing complex functions of position in the space, and 9fte denoting the real part. 
Then, the time-averaged Lorentz force over a time interval T large compared to 
2n/co [Born 99b] is 



<FW) = ^/_ 



l r T / 2 

dt 

T/2 



(p + p*) • V(E + E*) + ± + x (B + B*) 



(2) 



where * denotes complex conjugate. On substituting in Equation (2) E, B, and 
p by their time harmonic expressions given above and performing the integral, 
one obtains for each i'^-Cartesian component of the force 



WW> -^{ R 2gO + ^}. (3, 

In Equation (3) j = 1,2,3, is the completely antisymmetric Levy-Civita 
tensor. On using the Maxwell equation B = (c/iu)V x E and the relationships 
p = «E and dp/dt = —iup, a being the particle polarizability, Equation (3) 
transforms into 



1 f / dE*(r) f)F*\ s \ 

Since €jkiei mn = 8jm8kn—$jnfikm, one can finally express the time-averaged Lorentz 
force on the small particle as [Chaumet 00c] 



(^(r^^foE^}. (5) 



Equation (5) constitutes the expression of the time-averaged force on a particle 
in an arbitrary time-harmonic electromagnetic field. 

For a dipolar particle, the polarizability is [Draine 88] 



In Equation (6) «o is given by: ao = a 3 (e — l)/(e + 2), e = e2/eo being the 
dielectric permittivity contrast between the particle, 62, and the surrounding 
medium, eo; and k = ^/eoko, ^0 — oj/c. For ka <C 1, one can approximate a by: 
a = a (l + |iA; 3 |Q;o| 2 ). The imaginary part in this expression of a constitutes the 
radiation-reaction term. 

The light field can be expressed by its paraxial form, e.g., it is a beam or a 
plane wave, either propagating or evanescent, so that it has a main propagation 
direction along k, the light electric vector will then be described by 

E(r) = E (r)exp(ik-r). (7) 
Substituting Equation (7) into Equation (5), one obtains for the force 

(F) = ^e{a}V|E | 2 + ^m{a}\E \ 2 - ^m{a}%m{E • VE*}, (8) 

where 9m denotes imaginary part. The first term is the gradient force acting on 
the particle, whereas the second term represents the radiation pressure contri- 
bution to the scattering force that, on substituting the above approximation for 
a, namely, a = «o (1 + ^ik 3 \ «o| 2 ), can also be expressed for a Rayleigh particle 
(ka <C 1) as [van de Hulst 81] (|E | 2 /87r)Ck/A;, where C is the particle scattering 
cross section given by C = (8/3)7r/c 4 |a | 2 - Notice that the last term of Equa- 
tion (8) is only zero when either a or E is real. (This is the case for a plane 
propagating or evanescent wave but not for a beam, in general.) 

3 Force on a Dipolar Particle due to 
an Evanescent Wave 

Let the small particle be exposed to the electromagnetic field of an evanescent 
wave, whose electric vector is E = Texp(— qz) exp(iK-R), where we have written 
r = (R, z) and k = (K, k z ), K and k z satisfying K 2 + k 2 = k 2 , k 2 = u 2 e /c 2 , 

with k z = iq = i \J K 2 — k^ . This field is created under total internal reflection at 
a flat interface (z = constant, below the particle) between two media of dielectric 
permittivity ratio e /ei (see also inset of Figure 2(a)). The incident wave, either 
s or p polarized, (i.e., with the electric vector either perpendicular or in the plane 
of incidence: the plane formed by the incident wavevector kj at the interface and 
the surface normal z) enters from the denser medium at z < 0. The particle is in 
the medium at z > 0. Without loss of generality, we shall choose the incidence 
plane as OXZ, so that K = (K, 0). Let T± and Ty be the transmitted amplitudes 
into z > for s and p polarizations, respectively. The electric vector is: 
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E = (0, 1, 0)Tj_ exp(iKx) exp(-qz), 



(9) 



for s polarization, and 



TIi 

E = (-iq,0,K)-£exp(iKx)exp(-qz). (10) 

K 

for p polarization 

By introducing the above expressions for the electric vector E into Equa- 
tion (8), we readily obtain the average total force on the particle split into the 
scattering and gradient forces. The scattering force is contained in the OXY- 
plane (that is, the plane containing the propagation wavevector of the evanescent 
wave), namely, 



(Fx) = yi^S™{«} exp(-2qz); (11) 
For the gradient force, which is purely directed along OZ, one has 

(F z ) = -H^ e { a } exp(-2^). (12) 

In Equations (11) and (12) T stands for either Tj_ or Tj|, depending on whether 
the polarization is s or p, respectively. 

For an absorbing particle, on introducing Equation (6) for a into Equa- 
tions (11) and (12), one gets for the scattering force 



= H-k^-w rrw¥4r • (13) 



and for the gradient force 



^=-^ Twm^ eM - 2qz) - (14) 

It should be remarked that, except near resonances, in general ^sm{ao} is a 
positive quantity and therefore the scattering force in Equation (13) is positive 
in the propagation direction K of the evanescent wave, thus pushing the particle 
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parallel to the surface, whereas the gradient force Equation (14) is negative or 
positive along OZ, therefore, attracting or repelling the particle towards the 
surface, respectively, according to whether 9fte{a} > or 9fte{a} < 0. The 
magnitudes of these forces increase with the decrease of distance to the interface 
and it is larger for p polarization since in this case the dipoles induced by the 
electric vector at both the particle and the surface are oriented parallel to each 
other, thus resulting in a smaller interaction than when these dipoles are induced 
in the OXZ-plane (s-polarization) [Alonso 68]. 

In particular, if ka <C 1, Equation (13) becomes 



(F x ) 



\T\ 



-K exp(—2qz) 



3rv 

a <sm 



e- 1 



e + 2 



(15) 



The first term of Equation (15) is the radiation pressure of the evanescent wave on 
the particle due to absorption, whereas the second term corresponds to scattering. 
This expression can be further expressed as 



(F x } = ^jexp(-2qz)C ext . (16) 
where the particle extinction cross section C ext has been introduced as 



(17) 

Notice that Equation (17) coincides with the value obtained from Mie's theory 
for small particles in the low-order expansion of the size parameter ka of the 
extinction cross section [van de Hulst 81]. 

Although the above equations do not account either for multiple scattering 
as described by Mie's theory for larger particles or for multiple interactions of 
the wave between the particle and the dielectric surface, they are useful to un- 
derstand the fundametals of the force effects induced by a single evanescent wave 
on a particle. It should be remarked, however, that as shown at the end of this 
section, once Mie's theory becomes necessary, multiple scattering with the surface 
demands that its contribution be taken into account. 

Figure 1 shows the evolution of the scattering and gradient forces on three 
kinds of particles, namely, glass (e 2 = 2.25), silicon (e 2 = 15 + iO.lA) and gold 
(e 2 = —5.65 + i0.75), all of radius a = 60 nm, as functions of the gap distance 
d between the particle and the surface at which the evanescent wave is created. 
The illuminating evanescent wave is due to refraction of an incident plane wave 
of power P = ^0-\A\ 2 = 1.9 x 10~ 2 mW/jim 2 , equivalent to 150 mW over a 



C ext = Anka^m f 1 + — k 4 a 6 

U + 2J 3 e + 2 
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Figure 1: Forces in the Z direction and in the X direction (as insets) acting on 
a sphere with radius a = 60 nm, in the dipolar approximation. The angle of 
incidence is 9$ = 42°, larger than the critical angle 8 C = 41.8° (for the glass-air 
interface). A = 632.8 nm. Solid lines: s polarization, dashed lines: p polarization. 
The sphere material is: (a): glass, (b): silicon, (c): gold. 

circular section of radius 50 fim, on a glass-air interface at angle of incidence 
#o = 42° and A = 632.8 nm (the critical angle is 6 C = 41.8°), both at s and 
p polarization (electric vector perpendicular and parallel, respectively, to the 
incidence plane at the glass-air interface, namely, | Tj_ | 2 = 4ei cos 2 9 \A\ 2 /(e 1 — 1), 
m | 2 = 4 £l cos 2 O \A\ 2 /[( ei - 1)((1 + ei) sin 2 9 - 1)]). 

These values of forces are consistent with the magnitudes obtained on similar 
particles by applying Maxwell's stress tensor (to be discussed in Section 4.1) 
via Mie's scattering theory [Almaas 95]. However, as shown in the next section, 
as the size of the particle increases, the multiple interaction of the illuminating 
wave between the particle and the substrate cannot be neglected. Therefore, 
the above results, although of interpretative value, should be taken with care at 
distances smaller than 10 nm, since in that case multiple scattering makes the 
force stronger. This will be seen next. 



4 Influence of Interaction with the Substrate 

Among the several studies on forces of evanescent waves over particles, there are 
several models which calculate the forces from Maxwell's stress tensor on using 
Mie's theory to determine the scattered field in and around the sphere, with- 
out, however, taking into account the multiple scattering between the interface 
at which the evanescent wave is created and the sphere [Almaas 95, Chang 94, 
Walz 99]. We shall next see that, except at certain distances, this multiple inter- 
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action cannot be neglected. 

The equations satisfied by the electric and magnetic vectors in a non-magnetic 
medium are 



VxVxE-A; 2 E = 4vrA; 2 P, (18) 
V x V x H — k 2 H = -iAnkV x P, (19) 

where P is the polarization vector. The solutions to Equations (18) and (19) are 
written in integral form as 



E(r 
H(r 



) = k 2 j d 3 r' P(r') -g(ry), (20) 
= -ik J d 3 r V x P(r') • <?(r, r'), (21) 



In Equations (20) and (21) Q(r, r') is the outgoing Green's dyadic or field created 
at r by a point dipole at r'. It satisfies the equation 



V x V x g{r, r') - k 2 g(r, r') = 4vr5(r - r')l. (22) 
Let's introduce the electric displacement vector D = E + 4nP. Then, 



VxVxE = VxVxD - 4ttV x V x P. (23) 

Using the vectorial identity V x V x D = V(V • D) — V 2 D, and the fact that in 
the absence of free charges, V • D = 0, it is easy to obtain 



V x V x E = -V 2 D - 4tt [V(V • P) - V 2 P] (24) 
= -V 2 E-4ttV(V-P), (25) 
which straightforwardly transforms Equation (18) into: 

V 2 E + fc 2 E = -4vr[A; 2 P + V(V-P)], (26) 

whose solution is 



E(r) = k 2 J d 3 r' [P + V(V • P)] (r') G(r, r') . (27) 

In a homogenous infinite space the function G of Equation (27) is Go(r,r') = 
exp(ik\r — r'|)/|r — r'|, namely, a spherical wave, or scalar Green's function, 
corresponding to radiation from a point source at r'. 

To determine Q, we consider the case in which the radiation comes from a 
dipole of moment p, situated at r ; the polarization vector P is expressed as 

P(r) = p<J(r-r ). (28) 

Introducing Equation (28) into Equation (27) one obtains the well known expres- 
sion for the electric field radiated by a dipole 



E(r) = A; 2 pV(p • V) y ^ { _ lJ . (29) 
On the other hand, if Equation (28) is introduced into Equation (20) one obtains 

E(r) = £; 2 p-C?o(r,r'). (30) 

On comparing Equations (29) and (30), since both give identical value for E, we 

get 

k 2 P(v') ■ <? (r, r') = [A; 2 pV(p • V)]G (r, r'), (31) 
i.e., the tensor Green's function in a homogenous infinite space is 

So(r,r')= (z+lvv)Go(ry). (32) 

A remark is in order here. When applying Equation (32) in calculations one 
must take into account the singularity at r = r', this is accounted for by writing 

as [Yaghjian 80]: 

So(r, O = V [ (l + ^ VV) G (r, r')] - ^S(r - r')L v . (33) 

In Equation (33) V represents the principal value and L v is a dyadic that describes 
the singularity and corresponds to an exclusion volume around r = r', on whose 
shape it depends [Yaghjian 80]. 
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4.1 The Coupled Dipole Method 



Among the several methods of calculating multiple scattering between bodies 
of arbitrary shape (e.g., transition matrix, finite-difference time domain, integral 
procedures, discrete dipole approximation, etc.) we shall next address the coupled 
dipole method (Purcell and Pennipacker [Purcell73]). This procedure is specially 
suitable for multiple scattering between a sphere and a flat interface. 

Let us return to the problem of determining the interaction of the incident 
wave with the substrate and the sphere. The scattered electromagnetic field is 
obtained from the contribution of all polarizable elements of the system under 
the action of the illuminating wave. The electric vector above the interface is 
given by the sum of the incident field Ej and that expressed by Equation (20) 

with the dyadic Green function Q being given by 

£(r,r') = So(r,r') + S s (r,r'). (34) 

In Equation (34) Q is given by Equation (32) and, as such, it corresponds to the 

field created by a dipole in a homogeneous infinite space. On the other hand, Q s 
represents the field from the dipole after reflection at the interface. 

The polarization vector P is represented by the collection of N dipole moments 
Pj corresponding to the N polarizable elements of all materials included in the 
illuminated system, namely, 

N 

PW=E#-r»). (35) 

3 

The relationship between the k th dipole moment and the exciting electric field 
is, as before, given by pk = aifcE(rfc), with at expressed by Equation (6). Then, 
Equations (27), (34) and (35) yield 

N «• 

E(r,) = k 2 J2 afc[£o(r,, r fc ) + £,(r„ r fc )] • E(r fc ). (36) 
k 

The determination of Q s either above or below the flat interface is discussed 
next (one can find more details in Ref. [Agarwal 75a]). Let us summarize the 
derivation of its expression above the surface. The field E in the half-space 
z > 0, from a dipole situated in this region, is the sum of that from the dipole in 
free space and the field E r produced on reflection of the latter at the interface. 
Taking Equation (30) into account, this is therefore 
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(b): a=100 nm 
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Figure 2: Normalized force in the Z direction acting on a glass sphere on a glass- 
vacuum interface. The angle of incidence 9q = 42° is larger than the critical angle 
6 C = 41.8°. A = 632.8 nm. Thin lines: S polarization, Thick lines: P polarization, 
(a): a = 10 nm, full line: dipole approximation, dashed line: CDM-A, dotted 
line: CDM-B. The inset shows the scattering geometry, (b): a = 100 nm, 
full line: calculation with CDM-B, dashed line: static approximation. (From 
Ref. [Chaumet 00a]). 



E(r) = A; 2 p-^o(r,r , ) + E r (r), z>0. (37) 

Both the spherical wave Go an d E r are expanded into plane waves. The former 
is, according to Weyl's representation [Banos 66a, Nieto-Vesperinas 91], 

j fOO (ft l{ 

Go(r ' r,) = 27 loo kJK) GXP[ * (K ' (R " R0 + hlz - A)] ■ (38) 

On the other hand, E r is expanded as an angular spectrum of plane waves 
[Nieto-Vesperinas 91, Mandel 95] 

/oo 
d 2 K A r (K)exp[i(K-R + fc a z)]. (39) 
-oo 

Introducing Equation (38) into Equation (34) one obtains a plane wave expan- 
sion for Qq. This gives the plane wave components A^(K) of the first term of 
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Equation (37). Then, the plane wave components A r (K) of the second term of 
Equation (37) are given by 



A r (K) = r(K)A fc (K), 



(40) 



In Equation (40) r(K) is the Fresnel reflection coefficient corresponding to the 
polarization of A^. The result is therefore that Q, Equation (34), is 



g(r,r') = — / d 2 KS (K).g(K,z,^-S(K)exp[iK-(R-R')], 



47T 2 J-oo 

where [Agarwal 75a, Agarwal 75b, Keller 93] 



(41) 



j / k x ky 

SfK) = — \ -k y k x 
A ! OK 



(42) 



and the dyadic g has the elements [Greffet 97] 



9u = 

922 = 

933 = 
+ 

912 = 

913 = 

931 = 



2e k 2 



% % *M*kf\z + ,')] + exp(i*J»|* - z'\) 



2k 



(0) 



e k ( z + 



i.(o) , k .W 



exp[ik®(z + z')\ + exp(ifc(°>|2 - z'\) 



, (43) 
(44) 



iK 2 



1t$kz ^kg _€okz > + €\k 
1 



C>M - " ' expfifcf (z + z')\ - exp(ikW \z - z'\) 



(o) 



0, 



-iK 

iK 
2^kl 



5(z-z'), 



2 ^ "I - £\kz 

e k^ - Cl Af ) 



exp[i^ 0) (^ + *')] " exp(ikW\z - z'\) 



e kz + €ik z 

(0) 

— - -^eMik^iz + z'^+eMik^z-z'l) 

€ok z + €\k z 



(45) 
(46) 
(47) 

(48) 



We have used k[^ = iqj = (K 2 — e.,/^) 1 / 2 , j = ,1, and k = u/c. To determine 
the force acting on the particle, we also need the magnetic field. This is found by 
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the relationships B(r) = —i/kV x E(r). Then the time-averaged force obtained 
from Maxwell's stress tensor T [Stratton 41, Jackson 75] is 



(F> = \Sr (T(r)) 



■ n. 



(49) 



Equation (49) represents the flow of the time-average of Maxwell's stress tensor 

(T) across a surface S enclosing the particle, n being the local outward normal. 
The elements T a p are [Jackson 75] 



For dipolar particles, one can use instead of Maxwell's stress tensor the expression 
given by Equation (5) directly. In fact, for dielectric spheres of radii smaller 
than 5 x 10~ 2 A there is no appreciable difference between using Equation (5) or 
Equation (50) except at distances from the flat substrate smaller than 10~ 3 A. 

Figure 2 shows the normalized Z-force for two glass particles (e = 2.25) at 
A = 632.8 nm, one with a = 10 nm (Figure 2(a)), and other with a = 100 
nm (Figure 2(b)), the flat interface is illuminated from the dielectric side at 
#o = 42°, (the critical angle is 9 C = 41.8°). Two calculation procedures arc 
shown: a multiple scattering evaluation of the field via Equations (36)-(48) and 
then, either use of Equation (5), integrated over all induced dipoles (CDM-B), 
or Equation (49) (CDM-A). The normalization of the forces has been carried out 
by dividing them by exp(— 2qz). Thus, as seen in these curves, the force tends, 
as d increases, to the constant value given by Equation (12): — (\T\ 2 /2)q$te{a}. 
The incident power is 1.19 mW distributed on a surface of 10 fim 2 , then the 
force on a sphere of a = 10 nm is 2.7991 x 10~ 10 pN [Chaumet 00a]. We see, 
therefore, the effect on the vertical force of the multiple interaction of the scattered 
wave with the substrate: as the particle gets closer to the flat interface at which 
the evanescent wave is created, the magnitude of the attractive force increases 
beyond the value predicted by neglecting this interaction. As the distance to the 
surface grows, the force tends to its value given by Equation (12) in which no 
multiple scattering with the substrate takes place. Also, due to the standing wave 
patterns that appear in the field intensity distribution between the sphere and the 
substrate, the magnitude of this force oscillates as d varies. This is appreciated 
for larger particles (Figure 2(b)) except for very small particles (Figure 2(a)), 
whose scattering cross section is large enough to produce noticeable interferences. 
On the other hand, the horizontal force on the particle is of the form given by 
Equation (11) and always has the characteristics of a scattering force. 

As regards a metallic particle, we notice that 9fte{a} may have negative values 
near plasmon resonances (Figure 3(a), where we have plotted two models: that 



i r 



E a E; + B a B*p - I(E • E* + B • B*)5 a ^ , (a, = 1, 2, 3). 



8tt L 



(50) 
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Figure 3: (a): From top to bottom: the first three curves represent the polar- 
izability of a silver sphere with radius a = 10 nm versus the wavelength. The 
fourth curve is the force on this particle in free space. Plain line: Mie calculation, 
dashed line: polarizability of Eq. (6), symbol +: Dungey and Bohren's polar- 
izability [Dungey 91]. (b): Force along the Z direction on a silver sphere with 
a = 100 nm versus distance d with 6 = 50° for the following wavelengths: Plain 
line: A = 255 nm, dashed line: A = 300 nm, and dotted line: A = 340 nm. Thin 
lines: S polarization, thick lines: P polarization. (From Ref. [Chaumet 00b]). 

of Draine [Draine 88] and that of Ref. [Dungey 91] (see also [Chaumet 00c])) 
and thus the gradient force, or force along OZ, may now be repulsive, namely, 
positive (Figure 3(b)) [Chaumet 00c]. We also observe that this force is larger 
at the plasmon polariton resonance excitation (A = 350 nm). We shall later 
return to this fact (Section 5.1). We next illustrate how, no matter how small the 
particle is, the continuous approach to the surface makes the multiple scattering 
noticeable. 



4.2 Corrugated Surfaces: Integral Equations for 
Light Scattering from Arbitrary Bodies 

At corrugated interfaces, the phenomenon of TIR is weakened, and the contri- 
bution of propagating components to the transmitted field becomes important, 
either from conversion of evanescent waves into radiating waves, or due to the 
primary appearance of much propagating waves from scattering at the surface 
defects. The size of the asperities is important in this respect. However, TIR ef- 
fects are still strong in slightly rough interfaces, namely, those for which the defect 



14 



size is much smaller than the wavelength. Then the contribution of evanescent 
components is dominant on small particles (i.e., of radius no larger than 0.1 
wavelengths). On the other hand, the use of such particles as probes in near-field 
microscopy may allow high resolution of surface details as they scan above it. We 
shall next study the resulting force signal effects due to corrugation, as a model 
of photonic force microscopy under TIR conditions. 

When the surface in front of the sphere is corrugated, finding the Green's 
function components is not as straightforward as in the previous section. We 
shall instead employ an integral method that we summarize next. 

Let an electromagnetic field, with electric and magnetic vectors E^ nc ^(r) and 
jj( mc )(r), respectively, be incident on a medium of permittivity e occupying a 
volume V, constituted by two scattering volumes V\ and V 2 , each being limited by 
a surface Si and S 2 , respectively. Let r K be the position vector of a generic point 
inside the volume Vj, and by r > that of a generic point in the volume V , which 
is outside all volumes Vj. The electric and magnetic vectors of a monochromatic 
field satisfy, respectively, the wave equations, i.e., Equations (18) and (19). 

The vector form of Green's theorem for two vectors P and Q well behaved in 
a volume V surrounded by a surface S reads [Morse 53] 



f d 3 r (Q-VxVxP-P-VxVxQ) = 

d 2 r (P x V x Q - Q x V x P) • n, (51) 

with n being the unit outward normal. 

Let us now apply Equation (51) to the vectors P = Q(r, r') • C , (C being a 
constant vector) and Q = E(r). Taking Equations (18) and (22) into account, 
we obtain 



J d 3 r' E(r')<5(r - r') = k 2 J d 3 r' P(r') • £(r, r') - -^S e (r), (52) 



where S P is 



S e (r) = V x V x jf dV (E(r')^l^ - G(r,v')^^j . (53) 

Equation (53) adopts different forms depending on whether the points r and r' 
are considered in V or in V. By means of straightforward calculations one obtains 
the following: 
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• If r and r' belong to any of the volumes Vj, (j = 1,2), namely, V becomes 
either of the volumes Vf 

E(r<) = k 2 J d 3 r' P(r') • S^V) - ^Sf ) (r < ), (54) 

where 

S M( r <) = VxVx 

/ s/ v( Ej „ M ^- G(r <,^). w> 

In Equation (55) E in represents the limiting value of the electric vector on 
the surface Sj taken from inside the volume Vj. Equation (54) shows that 
the field inside each of the scattering volumes Vj does not depend on the 
sources generated in the other volumes. 

• If r belongs to any of the volumes Vj, namely, V becomes Vj, and r' belongs 
to V: 



= S ext (r<). (56) 

In Equation (56) S ext is 

SUr<) = £sfV)-Soo(r<), (57) 

j 

where 

s (-*)( r <) = VxVx 

(^'(^-^•^)- m 

In Equation (58) the surface values of the electric vector are taken from the 
volume V. The normal n now points towards the interior of each of the 
volumes Vj. 
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Also, Soo has the same meaning as Equation (58), the surface of integra- 
tion now being a large sphere whose radius will eventually tend to infinity. 
It is not difficult to see that — in Equation (57) is 4n times the inci- 
dent field 'E^ nc \r < ) (c/. Refs. [Nieto-Vesperinas 91] and [Pattanayak 76b, 
Pattanayak 76a]). Therefore Equation (56) finally becomes 



= E^ C V<) + E Sf ut \r<). (59) 

3 

Note that when Equation (59) is used as a non-local boundary condition, 
the unknown sources to be determined, given by the limiting values of 
E(r') and <9E(r')/<9n on each of the surfaces Sj, (cf. Equation (58)), appear 
coupled to those corresponding sources on the other surface Sk, k ^ j. 

Following similar arguments, one obtains: 

• For r belonging to V and r' belonging to either volume Vj, (j — 1, 2) namely, 
V becoming Vj 

= k 2 j d 3 r' P(r') • 5(r>, r) - i-Sf } (r>), (60) 

with Sj- m ^ given by Equation (55), this time evaluated at r > . 

• For both r and r' belonging to V 



E(r>) = E^ c )(r>) + i- ^ s5° u * ) (r > ), (61) 

j 

Hence, the exterior field is the sum of the fields emitted from each scattering 
surface Sj (j = 1, 2) with sources resulting from the coupling involved in 
Equation (60). 

One important case corresponds to a penetrable, optically homogeneous, isotropic, 
non-magnetic and spatially nondispersive medium (this applies for a real metal 
or a pure dielectric). In this case, Equations (54) and (60) become, respectively, 



E(r<) = --L-V x V x 
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= E (inc) (r < ) + V x V x 

?/ s /v( EM ^- G(r <,^), (as) 

whereas Equations (56) and (61) yield 



/ s /v( Ei „ M ^^- G - (r >,^), (04) 
E(r>) = E (mc) (r>) + x V x 

In Equations (62) and (64) "m" means that the limiting values on the surface 
are taken from inside the volume Vf note that this implies for both and E in 
that k = k yfe. 

The continuity conditions 

n x [E m (r<) - E(r>)] = 0, n x [H m (r<) - H(r>)] = , (66) 

and the use of Maxwell's equations lead to (cf. Ref. [Jackson 75], Section 1.5, or 
Ref. [Born 99c], Section 1.1): 



Ei, 



ires 



(-) 



£( r )U( + ) 



(67) 



dE m (r) 



dn 



res^ 



dE(r) 



dn 



res 



(+) 



(6£ 



where Sj + ^ and Sj~' denote the surface profile when approached from outside or 
inside the volume Vj, respectively. Equations (67) and (68) permit to find both 
E and d'E/dn from either the pair Equations (64) and (65), or, equivalently, 
from the pair Equations (62) and (63), as both r> and r< tend to a point in 
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Sj. Then the scattered field outside the medium is given by the second term of 
Equation (65). 

In the next section, we apply this theory to finding the near-field distribution 
of light scattered from a small particle in front of a corrugated dielectric surface 
when illumination is done from the dielectric half-space at angles of incidence 
larger than the critical angle. The non-local boundary conditions that we shall 
use are Equations (64) and (65). 

5 Photonic Force Microscopy of Surfaces with 
Defects 

The Photonic Force Microscope (PFM) is a technique in which one uses a probe 
particle trapped by a tweezer trap to image soft surfaces. The PFM [Ghislain 93, 
Florin 96, Wada 00] was conceived as a scanning probe device to measure ultra- 
small forces, in the range from a few to several hundredths pN/nm with laser pow- 
ers of some mW, between colloidal particles [Crocker 94], or in soft matter com- 
ponents such as cell membranes [Stout 97] and protein or other macromolecule 
bonds [Smith 96] . In such a system, a dielectric particle of a few hundred nanome- 
ters, held in an optical tweezer [Ashkin 86, Clapp 99, Sugiura 93, Dogariu 00], 
scans the object surface. The spring constant of the laser trap is three or four 
orders of magnitude smaller than that of AFM cantilevers, and the probe posi- 
tion can be measured with a resolution of a few nanometers within a range of 
some microseconds [Florin 96]. As in AFM, surface topography imaging can be 
realized with a PFM by transducing the optical force induced by the near field 
on the probe [Horber 01]. As in near-field scanning optical microscopy (NSOM 3 ) 
[Pohl 93], the resolution is given by the size of the particle and its proximity to 
the surface. It is well known, however [Nieto-Vesperinas 91, Greffet 97] that mul- 
tiple scattering effects and artifacts, often hinder NSOM images so that they do 
not bear resemblance to the actual topography . This has constituted one of the 
leading basic problems in NSOM [Hecht 96]. Numerical simulations [Lester 01, 
Arias-Gonzalez 01c, Arias- Gonzalez 99, Arias-Gonzalez 00, Arias-Gonzalez Ola] 
based on the theory of Section 4.2 show that detection of the optical force on the 
particle yields topographic images, and thus they provide a method of prediction 
and interpretation for monitoring the force signal variation with the topography, 
particle position and illumination conditions. This underlines the fundamentals 
of the PFM operation. An important feature is the signal enhancement effects 
arising from the excitation of Mie resonances of the particle, which we shall dis- 
cuss next. This allows to decrease its size down to the nanometric scale, thus 
increasing resolution both of force magnitudes and spatial details. 

3 NSOM is also called SNOM, abbreviation for scanning near-field optical microscopy. 



19 



5.1 Nanoparticle Resonances 

Electromagnetic eigenmodes of small particles are of importance in several areas 
of research. On the one hand, experiments on the linewidth of surface plasmons 
in metallic particles [Klar 88] and on the evolution of their near fields, both 
in isolated particles and in arrays [Krenn 99], seek a basic understanding and 
possible applications of their optical properties. 

Mie resonances of particles are often called morphology-dependent resonances 
(MDR). They depend on the particle shape, permittivity, and the size parameter: 
x = 2na/\. In dielectric particles, they are known as whispering-gallery modes 
(WGM) [Owen 81, Barber 82, Benincasa 87, Hill 88, Barber 88, Barber 90]. On 
the other hand, in metallic particles, they become surface plasmons (SPR), com- 
ing from electron plasma oscillations [Raether 88] . All these resonances are associ- 
ated to surface waves which exponentially decay away from the particle boundary. 

Morphology-dependent resonances in dielectric particles are interpreted as 
waves propagating around the object, confined by total internal reflection, re- 
turning in phase to the starting point. A Quality factor is also defined as Q = 2n 
(Stored energy) / (Energy lost per cycle) = uq/5u, where Uq is the resonace fre- 
quency and Sou the resonance full width. The first theoretical studies of MDR 
were performed by Gustav Mie, in his well-known scattering theory for spheres. 
The scattered field, both outside and in the particle, is decomposed into a sum of 
partial waves. Each partial wave is weighted by a coefficient whose poles explain 
the existence of peaks at the scattering cross section. These poles correspond to 
complex frequencies, but true resonances (i.e., the real values of frequency which 
produce a finite for the coefficient peaks) have a size parameter value close to 
the real part of the complex poles. The imaginary part of the complex fre- 
quency accounts for the width of the resonance peak. MDRs are classified by 
three integer numbers: one related to the partial wave {order number), another 
one which accounts for the several poles that can be present in the same co- 
efficient {mode number), and a third one accounting for the degeneration of a 
resonance (azimuthal mode number). In the first experimental check at optical 
frequencies, the variation of the radiation pressure (due to MDR) on highly trans- 
parent, low- vapor-pressure silicone oil drops (index 1.4 — 1.53) was measured by 
Ashkin [Ashkin 77]. The drops were levitated by optical techniques and the in- 
cident beam was focused at either the edge or the axis of the particles showing 
the creeping nature of the surface waves. 

It is important to note, as regards resonances, the enhanced directional scat- 
tering effects such as the Glory [Bryant 66, Fahlen 68, Khare 68] found in water 
droplets. The Glory theory accounts for the backscattering intensity enhance- 
ments found in water droplets. These enhancements are associated with rays 
grazing the surface of the droplet, involving hundreds of circumvolutions (sur- 
face effects). Axial rays (geometrical effects) also contribute. They have been 
observed in large particle sizes (x > 10 2 ) and no Glory effects have been found 
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for sizes in the range x ~ 1. These backscattering intensity enhancements cannot 
be associated to a unique partial wave, but to a superposition of several partial 
waves. 

5.2 Distributions of Forces on Dielectric Particles 

Over Corrugated Surfaces Illuminated Under TIR 

We now model a rough interface separating a dielectric of permittivity e\ = 
2.3104, similar to that of glass, from air. We have addressed (Figure 4, left) the 
profile consisting of two protrusions described by z = h[exp(— (x — X ) 2 /<7 2 ) + 
exp(— {x + Xq) 2 /<j 2 )] on a plane surface z — 0. (It should be noted that in 
actual experiments, the particle is immersed in water, which changes the particle's 
relative refractive index weakly. But the phenomena shown here will remain, 
with the interesting features now occurring at slightly different wavelengths.) 
Illumination, linearly polarized, is done from the dielectric side under TIR (critical 
angle 6 C = 41.14°) at 9q — 60° with a Gaussian beam of half-width at half- 
maximum W = 4000 nm at wavelength A (in air). For the sake of computing 
time and memory, the calculation is done in two dimensions (2D). This retains the 
main physical features of the full 3D configuration, as far as multiple interaction 
of the field with the surface and the probe is concerned [Lester 99]. The particle 
is then a cylinder of radius a, permittivity 62, and axis OY, whose center moves 
at constant height z — d + a. Maxwell's stress tensor is used to calculate the 
force on the particle resulting from the scattered near-field distribution created 
by multiple interaction of light between the surface and the particle. Since the 
configuration is 2D, the incident power and the force are expressed in mW/nm 
and in pN/nm, respectively, namely, as power and force magnitudes per unit 
length (in nm) in the transversal direction, i.e., that of the cylinder axis. We 
shall further discuss how these magnitudes are consistent with three-dimensional 
(3D) experiments. 

A silicon cylinder of radius a = 60 nm in front of a flat dielectric surface with 
the same value of e\ as considered here, has a Mie resonance excited by the trans- 
mitted evanescent wave at A = 638 nm [e^ = 14.99 + i0. 14) [Arias- Gonzalez 00]. 
Those eigenmodes are characterized by n — 0, 1 — 1, for p polarization, and n — 1, 
/ = 1, for s polarization. We consider the two protrusion interface (Figure 4, left). 
Insets of Figures 4(a) and 4(b), corresponding to a p-polarized incident beam, 
show the electromagnetic force on the particle as it scans horizontally above the 
flat surface with two protrusions with parameters = 63.8 nm, h = 127.6 nm 
and X = 191.4 nm, both at resonant A and out of resonance (A = 538 nm). The 
particle scans at d — 132.6 nm. Inset (a) shows the force along the OX axis. 
As seen, the force is positive and, at resonance, it has two remarkable maxima 
corresponding to the two protrusions, even though they appear slightly shifted 
due to surface propagation of the evanescent waves transmitted under TIR, which 



21 



XJX 




-800 -600 -400 -200 200 400 600 800 
X (nm) 



Figure 4: Left figure: Scattering geometry. Insets: Force curves on a silicon 
cylinder with a = 60 nm scanned at d — 132.6 nm. (a) Horizontal force, (b) 
Vertical force. Solid line: A = 638 nm (on resonance). Broken line: A = 538 nm 
(off resonance). Thin solid line in (b): \H/H \ 2 at z = d + a in absence of particle. 
Peak value: \H/H \ 2 = 0.07. Bottom figures: Spatial distribution \H/H \ 2 in this 
configuration. The cylinder center is placed at: (0, 192.6) nm (Fig. 4(a)), and at: 
(191.4,192.6) nm (Fig. 4(b)). The wavelength (A = 638 nm) excites the (n, I) 
Mie resonance. (From Ref. [Lester 01]). 

produce the Goos-Hanchen shift of the reflected beam. The vertical force on the 
particle, on the other hand, is negative, namely attractive, (cf. Inset (b)), and 
it has two narrow peaks at x just at the position of the protrusions. The signal 
being again remarkably stronger at resonant illumination. Similar force signal 
enhancements are observed for s-polarization. In this connection, it was recently 
found that this attractive force on such small dielectric particles monotonically 
increases as they approach a dielectric flat inteface [Chaumet 00a]. 

It should be remarked that, by contrast and as expected [Nieto-Vesperinas 91, 
Greffet 97], the near field-intensity distribution for the magnetic field H, normal- 
ized to the incident one H , has many more components and interference fringes 
than the force signal, and thus, the resemblance of its image with the interface 
topography is worse. This is shown in Inset (b) (thin solid line), where we have 
plotted this distribution in the absence of particle at z = d + a for the same 
illumination conditions and parameters as before. This is the one that ideally the 
particle scan should detect in NSOM. 

It is also interesting to investigate the near-field intensity distribution map. 
Figures 4(a) and 4(b) show this for the magnetic field H in p polarization for 
resonant illumination, A = 638 nm, at two different positions of the cylinder, 
which correspond to x = and 191.4 nm, respectively. We notice, first, the 
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Figure 5: Insets: Force on a silicon cylinder with a = 200 nm scanned at d — 442 
nm. (a) A = 919 nm (on resonance), (b): A = 759 nm (off resonance). Solid 
line: Vertical force. Broken line: Horizontal force. Figures: \E/E \ 2 in this 
configuration. The cylinder center is placed at: (0,642) nm (Fig. 5(a)), and at: 
(638,642) nm (Fig. 5(b)). The wavelength (A = 919 nm) excites the (n,l) Mie 
resonance. (From Ref. [Lester 01]). 

strong field concentration inside the particle, corresponding to the excitation of 
the (n — 0, 1 — 1) -eigenmode. When the particle is over one bump, the variation 
of the near field intensity is larger in the region closer to it, this being responsible 
for the stronger force signal on the particle at this position. Similar results are 
observed for s-polarized waves, in this case, the (n — 1, 1 — l)-eigenmode of the 
cylinder is excited and one can appreciate remarkable fringes along the whole 
interface due to interference of the surface wave, transmitted under TIR, with 
scattered waves from both the protrusions and the cylinder. 

An increase of the particle size yields stronger force signals at the expense of 
losing some resolution. Figures 5(a) and 5(b) show the near electric field inten- 
sity distribution for s polarization at two different positions of a cylinder with 
radius a = 200 nm, the parameters of the topography now being a = 212.7 nm, 
h = 425.3 nm and X = ±638 nm. The distance is d — 442 nm. The resonant 
wavelength is now A = 919 nm (e^ = 13.90+i0.07). Insets of Figures 5(a) and 5(b) 
illustrate the force distribution as the cylinder moves along OX. The force peaks, 
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when the resonant wavelength is considered, are positive, because now, the scat- 
tering force on this particle of larger scattering cross section is greater than the 
gradient force. They also appear shifted with respect to the protrusion positions, 
once again due to surface travelling waves under TIR. There are weaker peaks, 
or absence of them at non-resonant A. Similar results occur for a p-polarized 
beam at the resonant wavelength A = 759 nm (e 2 = 13.47 + iO.OA). For both po- 
larizations the (n — 3, 1 — 1) Mie eigenmode of the cylinder is now excited. The 
field distribution is well-localized inside the particle and it has the characteristic 
standing wave structure resulting from interference between counterpropagating 
whispering-gallery modes circumnavegating the cylinder surface. It is remarkable 
that this structure appears as produced by the excitation of propagating waves 
incident on the particle [Arias- Gonzalez 00], these being due to the coupling of 
the incident and the TIR surface waves with radiating components of the trans- 
mitted field, which are created from scattering with the interface protrusions. 
Although not shown here, we should remark that illumination at non-resonant 
wavelengths do not produce such a field concentration within the particle, then 
the field is extended throughout the space, with maxima attached to the flat por- 
tions of the interface (evanescent wave) and along certain directions departing 
from the protrusions (radiating waves from scattering at these surface defects). 

Evanescent components of the electromagnetic field and multiple scattering 
among several objects are often difficult to handle in an experiment. However, 
there are many physical situations that involve these phenomena. In this section 
we have seen that the use of the field inhomogeneity, combined with (and pro- 
duced by) morphology-dependent resonances and multiple scattering, permit to 
imaging a surface with defects. Whispering-gallery modes in dielectric particles, 
on the other hand, produce also evanescent fields on the particle surface which 
enhance the strength of the force signal. The next section is aimed to study 
metallic particles under the same situation as previously discussed, exciting now 
plasmon resonances on the objects. 

5.3 Distributions of Forces on Metallic Particles 

Over Corrugated Surfaces Illuminated Under TIR 

Dielectric particles suffer intensity gradient forces under light illumination due 
to radiation pressure, which permit one to hold and manipulate them by means 
of optical tweezers [Ashkin 77] in a variety of applications such as spectroscopy 
[Sasaki 91, Misawa 92, Misawa 91], phase transitions in polymers [Hotta 98], and 
light force microscopy of cells [Pralle 99, Pralle 98] and biomolecules [Smith 96]. 
Metallic particles, however, were initially reported to suffer repulsive electromag- 
netic scattering forces due to their higher cross sections [Ashkin 92], although 
later [Svoboda 94] it was shown that nanometric metallic particles (with diame- 
ters smaller than 50 nm) can be held in the focal region of a laser beam. Further, 
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Figure 6: \H/H Q \ 2 , P polarization, from a silver cylinder with a = 60 nm im- 
mersed in water, on a glass surface with defect parameters X = ±191.4 nm, 
h = 127.6 nm and a = 63.8 nm, at distance d = 132.6 nm. Gaussian beam inci- 
dence with W = 4000 nm. 6(a): A = 387 nm (on resonance), 9 Q = 0°. 6(b): 
A = 387 nm (on resonance), 9 a = 66°. 6(c): A = 316 nm (off resonance), 
9 Q = 66°. 6(d): A = 387 nm (on resonance), O = 66°. The cylinder center 
is placed at (0, 192.6) nm in 6(a), 6(b) and 6(c), and at (191.4, 192.6) nm in 6(d). 
(From Ref. [Arias- Gonzalez 01c]). 

it was demonstrated in an experiment [Sasaki 00] that metallic particles illumi- 
nated by an evanescent wave created under TIR at a substrate, experience a 
vertical attractive force towards the plate, while they are pushed horizontally in 
the direction of propagation of the evanescent wave along the surface. Forces in 
the fN range were measured. 

Plasmon resonances in metallic particles are not so efficiently excited as mor- 
phology-dependent resonances in non-absorbing high-refractive-index dielectric 
particles (e.g., see Refs. [Arias- Gonzalez 01b, Arias-Gonzalez 00]) under incident 
evanescent waves. The distance from the particle to the surface must be very 
small to avoid the evanescent wave decay normal to the propagation direction 
along the surface. In this section we address the same configuration as before 
using water as the immersion medium. The critical angle for the glass-water 
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interface is 9 C = 61.28°. A silver cylinder of radius a at distance d + a from the 
flat portion of the surface is now studied. 

In Figure 6 we plot the near-field intensity distribution \H/H \ 2 correspond- 
ing to the configuration of inset in Figure 4. A silver cylinder of radius a = 60 
nm scans at constant distance d = 162.6 nm above the interface. The system 
is illuminated by a p-polarized Gaussian beam (W = 4000 nm) at 9 = 0° and 
A = 387 nm (e 2 = —3.22 + i0.70). The surface protrusions are positioned at 
Xo = ±191.4 nm with height h = 127.6 nm and a = 63.8 nm. Figure 6(a) shows 
the aforementioned distribution when the particle is centered between the protru- 
sions. The plasmon resonance is excited as manifested by the field enhancement 
on the cylinder surface, which is higher in its lower portion. At this resonant 
wavelength, the main Mie coefficient contributor is n = 2, which can also be de- 
duced from the interference pattern formed along the particle surface: the number 
of lobes must be In along this surface [Owen 81]. Figure 6(b) shows the same 
situation but with 9 = 66°. The field intensity close to the particle is higher 
in Figure 6(a) because in Figure 6(b) the distance d is large enough to oblit- 
erate the resonance excitation due to the decay of the evanescent wave created 
by TIR [Arias-Gonzalez 01b]. However, the field intensity is markedly different 
from the one shown in Figure 6(c), in which the wavelength has been changed to 
A = 316 nm (e 2 = 0.78 + i 1.07) so that there is no particle resonance excitation at 
all. Figure 6(d) shows the same situation as in Figure 6(b) but at a different X- 
position of the particle. In Figure 6(c), the interference in the scattered near-field 
due to the presence of the particle is rather weak; the field distribution is now 
seen to be mainly concentrated at low z as an evanescent wave travelling along 
the interface, and this distribution does not substantially change as the particle 
moves over the surface at constant z. By contrast, in Figures 6(b) and 6(d) the 
intensity map is strongly perturbed by the presence of the particle. As we shall 
see, this is the main reason due to which optical force microscopy is possible at 
resonant conditions with such small metallic particles used as nanoprobes, and 
not so efficient at non-resonant wavelengths. In connection with these intensity 
maps (cf. Figures 6(b) and 6(d)), we should point out the interference pattern 
on the left side of the cylinder between the evanescent wave and the strongly re- 
flected waves from the particle, that in resonant conditions behaves as a strongly 
radiating antenna [Arias- Gonzalez 01b, Arias-Gonzalez 00, Krenn 99]. This can 
also be envisaged as due to the much larger scattering cross section of the parti- 
cle on resonance, hence reflecting backwards higher intensity and thus enhancing 
the interference with the evanescent incident field. The fringe spacing is A/2 (A 
being the corresponding wavelength in water). This is explained as follows: The 
interference pattern formed by the two evanescent waves travelling on the surface 
opposite to each other, with the same amplitude and no dephasing, is propor- 
tional to exp(— 2kz) cos 2 (nik sm9 x), with k = {n\ sin 2 9 — n 2 ,) 1 / 2 . The distance 
between maxima is Ax = X/(2n 1 sin# )- For the angles of incidence used in this 
work under TIR (9 = 66° and 72°), sin# 0.9, and taking into account the 
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refractive indices of water and glass, one can express this distance as Ax ps A/2no- 
The quantity Ax is similar to the fringe period below the particle in Figure 6(a), 
now attributted to the interference between two opposite travelling plane waves, 
namely, the one transmitted through the interface and the one reflected back 
from the particle. 




Figure 7: Force on a silver cylinder with a = 60 nm immersed in water, scanned 
at constant distance d = 132.6 nm on a glass surface with defect parameters 
X = ±191.4 nm, h = 127.6 nm and a = 63.8 nm along OX. The incident field is 
a p-polarized Gaussian beam with W = 4000 nm and 9 = 66°. 7(a): Horizontal 
force. 7(b): Vertical force. Solid curves: A = 387 nm (on resonance), broken 
curves: A = 316 nm (off resonance). Thin lines in 7(b) show \H/H \ 2 (in arbitrary 
units), averaged on the perimeter of the cylinder cross section, while it scans the 
surface. The actual magnitude of the intensity in the resonant case is almost seven 
times larger than in the non-resonant one. (From Ref. [Arias-Gonzalez 01c]). 

Figure 7 shows the variation of the Cartesian components of the electromag- 
netic force (F x , Fig. 7(a) and F z , Fig. 7(d)) on scanning the particle at constant 
distance d above the interface, at either plasmon resonance excitation (A = 387 
nm, solid lines), or off resonance (A = 316 nm, broken lines). The incident beam 
power (per unit length) on resonance is 3.9320 x 10~ 6 mW/nm, and 3.9327 x 10~ 6 
mW/nm at A = 316 nm. The incidence is done with ap-polarized Gaussian beam 
of W = 4000 nm at 9 = 66°. It is seen from these curves that the force distribu- 
tions resembles the surface topography on resonant conditions with a signal which 
is remarkably larger than off-resonance. This feature is specially manifested in 
the Z component of the force, in which the two protrusions are clearly distin- 
guished from the rest of interference ripples, as explained above. Figure 7(b) also 
shows (thin lines) the scanning that conventional near field microscopy would 
measure in this configuration, namely, the normalized magnetic near field inten- 
sity, averaged on the cylinder cross section. These intensity curves are shown in 
arbitrary units, and in fact the curve corresponding to plasmon resonant condi- 
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Figure 8: \H/H \ 2 for P polarization for a silver cylinder with a = 200 nm im- 
mersed in water, on a glass surface with parameters X = ±638 nm, h = 425.3 
nm and a = 212.7 nm, at distance d = 442 nm. Gaussian beam incidence with 
W = 4000 nm. 8(a): A = 441 nm (on resonance), 9 Q = 0° and the cylinder 
center placed at (—1276,642) nm. 8(b): A = 441 nm (on resonance), 9 Q = 66° 
and the cylinder center placed at (1276,642) nm. 8(c): A = 316 nm (off res- 
onance), 9 Q = 66° and the cylinder center placed at (1276,642) nm. (From 
Ref. [Arias- Gonzalez 01c]). 



tions is almost seven times larger than the one off-resonance. The force curves 
show, on the one hand, that resonant conditions also enhance the contrast of 
the surface topography image. Thus, the images obtained from the electromag- 
netic force follows more faithfully the topography than that from the near field 
intensity. This is a fact also observed with other profiles, including surface-relief 
gratings. When parameter h is inverted, namely, the interface profile on the left in 
Fig. 4, then the vertical component of the force distribution presents inverted the 
contrast. On the whole, one observes from these results that both the positions 
and sign of the defect height can be distinguished by the optical force scanning. 

Figure 8 displays near-field intensity maps for a larger particle (a = 200 nm). 
Figure 8(a) corresponds to 9 = 0° and a resonant wavelength A = 441 nm 
(e2 = —5.65 + i0.75), with the particle being placed on the left of both protru- 
sions. Figure 8(b) corresponds to 9 = 66° (TIR illumination conditions), at the 
same resonant wavelength, the particle now being on the right of the protrusions. 
Figure 8(c) corresponds to 9 = 66° (TIR incidence), at the no resonant wave- 
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Figure 9: Force on a silver cylinder with a = 200 nm immersed in water, 
scanned at constant distance on a glass surface with parameters X = ±638 
nm, h = 425.3 and a = 212.7 nm along OX. The incident field is a p- 
polarized Gaussian beam with W = 4000 nm and 9 = 66°. 9(a): Horizontal 
force. 9(b): Vertical force. Solid curves: A = 441 nm (on resonance), broken 
curves: A = 316 nm (off resonance). Thin solid curves: A = 441 nm (on reso- 
nance) at 6*o = 72°, thin broken curves: A = 316 nm (off resonance) at 9 = 72°. 
(From Ref. [Arias- Gonzalez 01c]). 

length A = 316 nm (e 2 = 0.78 + «1.07) , the particle being placed at the right 
of the protrusions. The incident beam is p-polarized with W = 4000 nm. The 
surface protrusions are positioned at X Q = ±638 nm with height h = 425.3 nm 
and a = 212.7 nm. All the relevant size parameters are now comparable to the 
wavelength, and hence to the decay length of the evanescent wave. That is why 
now the plasmon resonance cannot be highly excited. When no resonant wave- 
length is used, the intensity interference fringes due to the presence of the particle 
are weaker. On the other hand, Figure 8(a) shows the structure of the near-field 
scattered under 6 = 0°. There are three objects that scatter the field: the two 
protrusions and the particle. They create an inteference pattern with period A/2 
(with A being the wavelength in water). Besides, the particle shows an intefer- 
ence pattern around its surface due to the two counterpropagating plasmon waves 
which circumnavegate it [Arias- Gonzalez 01b, Arias- Gonzalez 00]. The number 
of lobes along the surface is nine, which reflects that the contribution to the field 
enhancement at this resonant wavelength comes from Mie's coefficients n = 5 
and n — 4. Figure 8(b) shows weaker excitation of the same plasmon resonance 
under TIR conditions. Now, the interference pattern at the incident side of the 
configuration is also evident. This pattern again has a period A/2 (A being the 
wavelength in water). If non-resonant illumination conditions are used, the parti- 
cle is too far from the surface to substantially perturb the transmitted evanescent 
field, then the intensity distribution of this field remains closely attached to the 
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interface, and it is scattered by the surface protrusions. The field felt by the 
particle in this situation is not sufficient to yield a well-resolved image of the 
surface topography as shown next for this same configuration. 

Figure 9 shows the components of the force (F x , Figure 9(a) and F z , Fig- 
ure 9(b)) for either plasmon excitation conditions (A = 441 nm, solid lines), 
or off-resonance (A = 316 nm, broken lines), as the cylinder scans at constant 
distance d above the surface. The incidence is done with a p-polarized Gaus- 
sian beam of W = 4000 nm at either 9q — 66° (thick curves) or 9q — 72° (thin 
curves). The incident beam power (per unit length) is 3.9313 x 10~ 6 mW/nm 
on resonance and 3.9327 x 10~ 6 mW/nm at A = 316 nm when 9 = 66°, and 
3.9290 x 10~ 6 mW/nm on resonance and 3.9315 x 10~ 6 mW/nm at A = 316 nm 
when 6 = 72°. As before, resonant conditions provide a better image of the 
surface topography making the two protrusions distinguishable with a contrast 
higher than the one obtained without plasmon excitation. The surface image 
corresponding to the force distribution is better when the protrusions (not shown 
here) are inverted because then the particle can be kept closer to the interface. 
Again, the curve contrast yielded by protrusions and grooves is inverted from 
each other. The positions of the force distribution peaks corresponding to the 
protrusions now appear appreciably shifted with respect to the actual protru- 
sions' position. This shift is explained as due to the Goos-Hanchen effect of the 
evanescent wave [Lester 01]. We observe that the distance between these peaks 
in the F z curve is aproximately 2X . This shift is more noticeable in the force 
distribution as the probe size increases 4 . Again, the F z force distribution has 
a higher contrast at the (shifted) position of the protrusions. The force signal 
with these bigger particles is larger, but the probe has to be placed farther from 
the surface at constant height scanning. This affects the strength of the signal. 
Finally, it is important to state that the angle of incidence (supposed to be larger 
than the critical angle C ) influences both the contrast and the strength of the 
force: the contrast decreases as the angle of incidence increases. At the same 
time, the strength of the force signal also diminishes. 

As seen in the force figures for both sizes of particles, most curves contain 
tiny ripples. They are due to the field intensity interference pattern as shown in 
Figures 6 and 8, and discussed above. As the particle moves, the force on it is 
affected by this interference. As a matter of fact, it can be noted in the force 
curves that these tiny ripples are mainly present at the left side of the particle, 
which is the region where stronger interference takes place. 

It is worth remarking, however, that these oscillations are less marked in 
the force distribution (c/. their tiny ripples), than in the near field intensity 
distribution, where the interference patterns present much higher contrast. 

As stated in the previous section, evanescent fields and multiple scattering are 
fruitful to extracting information from a detection setup. The latter is, at the 

4 For a better picture of this shift, see the grating case in Ref. [Lester 01]. 
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same time, somewhat troublesome as it cannot be neglected at will. This incove- 
nient is well-known in NSOM, but it is diminished in PFM, as remarked before. 
The smoother signal provided by the force is underlined by two facts: one is the 
averaging process on the particle surface, quantitatively interpreted from the field 
surface integration involved in Maxwell's stress tensor. Other is the local charac- 
ter of the force acting at each point of the particle surface. Metallic particles are 
better candidates as probes of PFM in comparison to dielectric particles, since 
the force signal is not only enhanced at resonance conditions, but it is also bigger 
and presents better resolution. However, dielectric particles are preferred when 
the distance to the interface is large, since then the weak evanescent field present 
at these distances presents better coupling to the whispering-gallery modes than 
to the plasmon surface waves of metallic particles. 

5.4 On the Attractive and Repulsive Nature of 

Vertical Forces and Their Orders of Magnitude 

The horizontal forces acting on the particle are scattering forces due to radiation 
pressure of both the incident evanescent wave and the field scattered by the 
protrusions, thus the forces are positive in all the cases studied. As for the vertical 
forces, two effects compete in determining their sign. First, is the influence of the 
polarizability [Chaumet 00b, Chaumet 00a], which depends on the polarization 
of the illumination. On the other hand, it is well known that an evanescent wave 
produces only gradient forces in the vertical direction. For silver cylinders, the 
force at wavelength A = 387 (e 2 = —3.22 + i0.70) and at A = 441 nm (e 2 = 
—5.65 + i0.75) must be attractive, while at A = 316 nm (e 2 = 0.78 + «1.07), the 
real part of the polarizability changes its sign, and so does the gradient force, thus 
becoming repulsive (on cylinders of not very large sizes, as here). However, in the 
cases studied here, not only the multiple scattering of light between the cylinder 
and the flat portion of the interface, but also the surface defects, produce scattered 
waves both propagating (into z > 0) and evanescent under TIR conditions. Thus, 
the scattering forces also contribute to the ^-component of the force. This affects 
the sign of the forces, but it is more significant as the size of the objects increases. 
In larger cylinders and defects (c/. Figure 9), the gradient force is weaker than 
the scattering force thus making F z to become repulsive on scanning at A = 
441 nm (plasmon excited). On the other hand, for the smaller silver cylinders 
studied (cf. Figure 7), the gradient force is greater than the scattering force at 
A = 387 nm (plasmon excited), and thus the force is attractive in this scanning. 
Also, as the distance between the particle and the surface decreases, the gradient 
force becomes more attractive [Chaumet 00b, Chaumet 00a]. This explains the 
dips and change of contrast in the vertical force distribution on scanning both 
protrusions and grooves. At A = 316 nm (no plasmon excited), both scattering 
and gradient forces act cooperatively in the vertical direction making the force 
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repulsive, no matter the size of the cylinder. For the silicon cylinder, as shown, the 
vertical forces acting under TIR conditions are attractive in absence of surface 
interaction (for both polarizations and the wavelengths used). However, this 
interaction is able to turn into repulsive the vertical force for S polarization at 
A = 538 nm, due to the scattering force. 

This study also reveals the dependence of the attractive or repulsive nature 
of the forces on the size of the objects (probe and defects of the surface), apart 
from the polarizability of the probe and the distance to the interface, when illu- 
mination under total internal reflection is considered. The competition between 
the strength of the scattering and the gradient force determines this nature. 

The order of magnitude of the forces obtained in the preceding 2D calcu- 
lations is consistent with that of forces in experiments and 3D calculations of 
Refs. [Pohl 93, Guntherodt 95, Depasse 92, Sugiura 93, Kawata 92, Dereux 94, 
Girard 94, Almaas 95, Novotny 97, Hecht 96, Okamoto 99, Chaumet 00b] and 
[Chaumet 00a]. Suppose a truncated cylinder with axial length L = 10 fim, and 
a Gaussian beam with 2W ~ 10 fim. Then, a rectangular section of L x 2W = 10 2 
fim 2 is illuminated on the interface. For an incident power P ~ 1 mW, spread 
over this rectangular section, the incident intensity is I ~ 10 -2 mW/fim 2 , and 
the force range from our calculations is F ~ 10~ 2 — 10 _1 pN. Thus, the forces ob- 
tained in Figures 9(b) and 9(d) are consistent with those presented, for example, 
in Ref. [Kawata 92]. 

6 Concluding Remarks & Future Prospects 

The forces exerted by both propagating and evanescent fields on small particles 
are the basis to understand estructural characteristics of time-harmonic fields. 
The simplest evanescent field that can be built is the one that we have illustrated 
on transmission at a dielectric interface when TIR conditions occur. Both di- 
electric and metallic particles are pushed along the direction of propagation of 
the evanescent field, independently of the size (scattering and absorption forces). 
By contrast, forces behave differently along the decay direction (gradient forces) 
on either dielectric and metallic particles, as studied for dipolar-sized particles, 
Section 2. The analysis done in presence of a rough interface, and with particles 
able to interact with it show that scattering, absorption and gradient forces act 
both in the amplitude and phase directions, when multiple scattering takes place. 
Moreover, the excitation of particle resonances enhances this interaction, and, at 
the same time, generates evanescent fields (surface waves) on their surface, which 
makes even more complex this mixing among force components. Thus, an analy- 
sis based on a small particle isolated is not feasible due to the high inhomogeneity 
of the field. 

It is, however, this inhomogeneity what provides a way to imaging a sur- 
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face with structural features, such as topography. The possibility offered by 
the combination of evanescent fields and Mie resonances is however not unique. 
As inhomgeneous fields (which can be analitically decomposed into propagat- 
ing and evanescent fields [Nieto-Vesperinas 91]) play an important role in the 
mechanical action of the electromagnetic wave on dielectric particles (either on 
or out of resonance), they can be used to operate at the nanometric scale on 
such entities, to assist the formation of ordered particle structures as for ex- 
ample [Burns 89, Burns 90, Antonoyiannakis 97, Antonoyiannakis 97, Malley 98, 
Bayer98, Barnes 02], with help of these resonances. Forces created by evanes- 
cent fields on particles and morphology-dependent resonances are the keys to 
control the optical binding and the formation of photonic molecules. Also, when 
a particle is used as a nanodetector, these forces are the signal in a scheme of 
photonic force microscopy as modeled in this article. It has been shown that the 
evanescent field forces and plasmon resonance excitations permit to manipulate 
metallic particles [Novotny 97, Chaumet 01, Chaumet 02], as well as to make 
such microscopy [Arias-Gonzalez 01c]. Nevertheless, controlled experiments on 
force magnitudes, both due to evanescent and propagating waves, are yet scarce 
and thus desirable to be fostered. 

The concepts released in this article open an ample window to investigate 
on soft matter components like in cells and molecules in biology. Most folding 
processes require small forces to detect and control in order not to alter them and 
be capable of actuating and extracting information from them. 

General Annotated References 

Complementary information and sources for some of the contents treated in this 
report can be found in next bibliography: 

• Electromagnetic optics: there are many books where to find the basis 
of the electromagnetic theory and optics. We cite here the most com- 
mon: [Jackson 75, Born 99c]. The Maxwell's Stress Tensor is analysed 
in [Jackson 75, Stratton 41]. The mathematical level of these textbooks 
is similar to the one in this report. 

• Mie theory can be found in [van de Hulst 81, Kerker 69, Bohren 83]. In 
these textbooks, Optics of particles is develop with little mathematics and 
all of them are comparable in contents. 

• Resonances can be understood from the textbooks before, but a more de- 
tailed information, with applications and the implications in many top- 
ics can be found in the following references: [Arias-Gonzalez 02c, Hill 88, 
Barber 90]. The last two references are preferently centred in dielectric 
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particles. The first one compiles some of the information in these two refer- 
ences and some other from scientific papers. Surface plasmons, on the other 
hand, are studied in depth in [Raether 88]. They are easy to understand 
from general physics. 

• Integral equations in scattering theory and angular spectrum representation 
(for the decomposition of time-harmonic fields in propagating and evanes- 
cent components) are treated in [Nieto-Vesperinas 91]. The mathematical 
level is similar to the one in this report. 

• The Coupled Dipole Method can be found in the scientific papers cited in 
Section 4.1. A More didactical reference is [Chaumet 98, Rahmani 00]. The 
information shed in this report on the CDM is extended in these references. 

• The dipolar approximation, in the context of optical forces and evanescent 
fields, can be complemented in scientific papers: [Gordon 73, Chaumet 00c, 
Arias-Gonzalez 02b] and in the monograph: [Novotny 00]. 

• A more detailed discussion on the sign of optical forces for dipolar par- 
ticles, as well as larger elongated particles, can be found in the scientific 
papers [Arias- Gonzalez 02a, Arias- Gonzalez 02b, Chaumet 00b]. 

• Monographs on NSOM and tweezers: [Nieto-Vesperinas 96, Sheetz 97]. 
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